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LIGHT  SCATTERING  BY  ROUGH  SURFACES  FULL  WAVE  SOLUTIONS 
1 .  Introduction 

1.1  Objectives,  Motivation  and  Background  for  this  Work 

Scattering  of  electromagnetic  vaves  in  media  consisting  of  random 
distributions  of  particles  such  as  atmospheric  aerosols,  smoke  and  dust 
has  been  investigated  extensively  using  the  equation  of  transfer 
(Chandrasekhar  1950,  Ishimaru  1978).  The  main  difficulties  in  setting 
up  the  equation  of  transfer  lie  in  the  determination  of  the  elements  of 
the  loci*  phase  (scattering)  matrix  and  the  extinction  cross  sections  for 
the  individual  particles.  Most  of  the  work  has  been  done  for  particles 
of  idealized  shapes  such  as  spheres  and  cylinders  (Ruck  et  al.  1970). 

In  many  physical  problems  of  interest  however,  the  individual 
scatterers  are  of  irregular  shapes  such  as  flakes,  spheres  and  cylinders 
with  random  rough  surfaces  (Scheurman  1980,  Bahar  and  Fitzvater  1983, 

Eahar  and  Chakrabarti  1985). 

In  this  work  the  results  derived  for  particles  with  idealized  shapes 
(spheres)  are  modified  to  account  for  the  random  surface  roughness  of  the 
particles.  To  this  end  the  full  wave  approach  (Bahar  1981)  is  used  to 
determine  the  rough  surface  contributions  to  the  like  and  cross  polarized 
scattering  cross  sections  (Bahar  and  Chakrabarti  1985),  and  the  elements 
of  the  phase  (scattering)  matrix  are  given  in  terms  of  a  weighted  sum  of 
the  Mie  solutions  and  the  diffuse  scattering  terms  due  to  the  particle 
surface  roughness.  The  extinction  cross  sections  are  determined  through  a 
judicious  use  of  the  forward  scattering  theorem  (Born  and  Wolf  1961+) 
and  the  very  perceptive  observation  that  for  large  scatterers  compared 
to  a  wavelength  the  forward  scattered  "shadow  forming  wave  is  the  same 
for  all  surfaces  which  have  the  same  shadowline"  (Morse  and  Feshbach  1953). 


Linearly  polarized  waves  are  assumed  to  be  incident  upon  a  parallel 
layer  consisting  of  a  random  distribution  of  irregular  shaped  particles. 


Particles  with  different  complex  dielectric  coefficients  and  different 
surface  roughness  are  considered.  The  rough  surface  height  is  characterized 
by  its  surface  height  spectral  density  functions  (the  Fourier  transforms 
of  the  non-Gaussian  surface  height  autocorrelation  functions).  Thus 
particles  with  a  large  range  of  roughness  scales  are  considered.  Special 
attention  is  given  to  determining  the  effects  of  varying  the  mean  square 
height  and  slope  of  the  particles. 

Beth  single  scatter  (first  order)  and  multiple  scatter  results  are 
presented  for  smooth  particles  and  for  particles  with  rough  surfaces. 

Layers  with  different  optical  thicknesses  are  considered  and  the  results 
are  presented  for  the  forward  and  backward  scattered  incoherent  diffuse 
scattering  intensities  as  well  as  the  co-polarized  and  cross  polarized 
intensities.  The  degree  of  polarization  of  the  reflected  and  transmitted 
intensities  are  also  determined.  The  Gaussian  quadrature  technique  is  used 
to  discretize  the  equation  of  transfer  and  the  matrix  characteristic 
value  technique  is  used  to  account  for  multiple  scattering  (ishimaru  and 
Cheung  1980).  This  work  is  particularly  relevant  to  problems  of  obscuration 
by  aerosols. 

1.2  Principal  Elements  of  the  Full  Wave  Approach 

As  indicated  in  Section  1.1,  the  full  wave  solutions  for  the  scattered 
like  and  cross  polarized  intensities  are  determined  using  the  full  wave 
approach  which  accounts  for  the  specular  point  (physical  optics)  scattering 
and  diffuse  scattering  in  a  self-consistent  manner.  Unlike  ihe  small 
perturbation  method  (Rice  1951,  Ruck  et  al.  1970,  Kiehl  et  al.  1980), 
the  full  wave  approach  is  not  restricted  to  surfaces  with  small  roughness 


parameters  B  =  l*k  <h  >  (vhere  <h  >  is  the  mean  square  height  and  kQ  is  the 
free  space  wavenumber  of  the  incident  electromagnetic  wave).  The  principal 
elements  of  the  full  wave  approach  (Bahar  1982)  are: 

(a)  The  electromagnetic  fields  are  expressed  in  terms  of  complete 
expansions  of  vertically  and  horizontally  polarized  waves.  These  include 
the  radiation  fields,  the  lateral  waves  and  the  surface  waves. 

(b)  Exact  boundary  conditions  are  imposed  at  the  irregular  surface. 

Cc )  Using  the  orthogonal  properties  of  the  basis  functions  appearing 

in  the  complete  expansions  of  the  fields.  Maxwell's  equations  are  integrated 
over  the  transverse  plane.  Green's  theorems  are  used  to  avoid  term-by-term 
differentiation  of  the  field  expansions. 

(d)  Maxwell's  equations  for  the  electromagnetic  fields  are  converted 
into  coupled  first  order  ordinary  differential  equations  for  the  forward 
and  backward  traveling  wave  amplitudes  which  are  only  functions  of  the 
variable  x.  (in  view  of  the  integration  in  the  transverse  plane  (y,z), 
the  telegraphists'  equations  are  only  fractions  of  x).  The  coupled 
equations  for  the  wave  amplitudes  are  referred  to  the  generalized 
telegraphists'  equations. 

(e)  A  variable  coordinate  system  that  conforms  with  the  local  features 
of  the  irregular  boundary  is  introduced  and  the  resulting  solutions  for 

the  scattered  fields  are  shown  to  be  invariant  to  coordinate  transformations 

(f)  Closed  form  second  order  iterative  solutions  for  the  radiation 
fields  are  obtained  from  the  telegraphists'  equations  on  neglecting 
multiple  scattering  from  the  rough  surface.  These  second  order  iterative 
solutions  account  for  wave  scattering  in  arbitrary  directions. 

(g)  The  full  wave  solutions  are  compared  with  earlier  geometric  optics 
physical  optics  and  perturbation  solutions.  The  suitability  of  the  two- 
scale  model  is  shown  to  be  limited. 


(i)  Evaluate  the  like  and  cross  polarized  scattering  cross  sections 
for  particles  of  irregular  shape.  The  particles  are  characterized 
"by  their  complex  permittivities  (unperturbed)  sizes,  rough  surface 
height  spectral  density  functions — (or  autocorrelation  functions). 
The  full  vave  approach  (see  Section  1.2),  is  used  to  account  for 
specular  point  and  diffuse  scattering  in  a  self-consistent  manner, 

(ii)  Evaluate  the  extinction  cross  sections  and  the  albedos  for  the 
irregular  particles  of  finite  conductivity. 

(iii)  Solve  the  transfer  equation  for  the  incoherent  specific  intensities 
(modified  Stokes  parameters),  and  the  co-polarized  and  cross 
polarized  intensities  for  linearly  polarized  vaves  at  optical 
and  infrared  frequencies  incident  upon  a  layer  of  randomly 
distributed  particles.  Compare  the  results  for  smooth  and  rough 
particles.  Present  both  single  scatter  and  multiple  scatter 
results . 

(iv)  Evaluate  the  degree  of  polarization  of  the  reflected  and  trans¬ 
mitted  intensities. 

(v)  Document  computer  programs  for  inclusion  in  the  final  renort . 


3.  Formulation  of  the  Problem 

The  transfer  equation  for  the  diffuse  specific  intensities  (Stokes 
parameters)  can  he  vritten  as  follows  in  matrix  notation  (Chandrasekhar 


1950,  Ishimaru  1978) 


V  ■  -[I]  +  /[S][l']dvi’d<f>'  +  [I.]  (1) 

in  which  T  is  the  optical  distance  in  the  z  direction  (normal  to  the 
plane  of  a  horizontal  layer  consisting  of  a  random  distribution  of 
particles  (see  Fig.  l) 

T  =  2p[ot]  =  Z  /  0  h(D)dD  (2) 

where  D  is  the  diameter  of  the  unperturbed  spherical  particle,  n(D) 
is  the  particle  size  distribution  and  o  is  the  extinction  coefficient 
(total  cross  section).  The  symbol  z[*]  denotes  integration  over  the  size 
distribution.  Since  the  unperturbed  particle  is  assumed  here  to  be  a 
sphere  of  diameter  D  and  the  surface  roughness  is  assumed  to  be  homogeneous 
and  isotropic,  diffuse  scattering  vanishes  in  the  forward  direction 
(Bahar  and  Chakrabarti  1985)  and  the  extinction  matrix  for  the  rough 
sphere  can  be  represented  by  a  scalar  quantity  (ishimaru  and  Cheung  1980). 
The  (lixl)  column  matrices  [i]  and  [i*]  are  the  incoherent  diffuse  specific 
intensity  matrices  (Stokes  parameters)  for  waves  scattered  by  particles 
in  the  direction  6  =  cos  and  $  and  for  waves  incident  in  the  direction 


0'  =  cos  y'  and  <p’  respectively.  The  elements  of  the  matrix  [i]  are 


the  modified  Stokes  parameters 


<E1E*> 

<EgE*> 

2Re<E]E*> 


2Im<E1E«> 


(3) 


in  which  the  symbol  <•>  denotes  the  statistical  average  and  the  symbol  * 
denotes  the  complex  conjugate.  An  exp(iut)  suppressed  time  dependent 


excitation  is  assumed.  The  vertically  and  horizontally  polarized 
components  of  the  electric  field  are  E^  and  E 2  respectively.  The  (hxk) 
scattering  (phase)  matrix  S  in  the  fixed  reference  coordinate  system 
(see  Fig.  l)  is  expressed  in  terms  of  the  scattering  matrix  S'  in  the 
scattering  plane  (vhich  contains  the  vave  normals  n^  and  n^  in  the  direction 
of  the  incident  and  scattered  waves)  through  the  following  transformation 
[S]  =  &f(-7r  +a)][S'][^a’)]  (k) 

in  which  [S']  is  the  weighted  sum  of  two  matrices 

[S']  =  |x(v-Sr>!2[sMle3  *  [SD]  (5) 

In  (k)  [of]  is  a  *txk  transformation  matrix  that  represents  the  specific 

intensities  relative  to  the  reference  plane  of  incidence  to  the  intensities 

relative  to  the  scattering  plane  and  back  from  the  scattering  plane  to  the 

reference  plane  of  scatter  (Chandrasekhar  1950,  Ishimaru  1978).  The 

matrix  [S...  ]  in  (5)  is  the  olase  matrix  for  the  unperturbed  spheres 
I-.ie 

based  cn  the  Kie  solution  (Ruck  et  al.  1970,  Ishimaru  1978).  Furthermore, 
in  (5)  x(”'*av.)  is  the  particle  random  rough  surface  height  characteristic 
function 

X(v*ar)  =  exp<ivarh>  (6) 

in  which  h  is  the  rough  surface  height  measured  in  the  direction  a^,  the  unit 
vector  in  the  radial  direction  (normal  to  the  surface  of  the  unperturbed 
sphere)  and 

v  =  k  (n^-n^ )  (7) 

o 

In  this  work  the  rough  surface  height  nrobability  density  function  p(h) 
and  joint  probability  function  p(h,h’)  are  assumed  to  be  Gaussian,  thus 

|x]2  =  exp(-v*ar<h2>)  (8) 

and  the  rough  surface  height  Joint  characteristic  function  is 


(9) 


=  <exp[ivar (h-h '  )]>  =  exp(v*ar<hh'>)  Jx| 

Thus  in  (5)  the  coefficient  ] x I  accounts  for  the  degradation  of  the 
scattering  by  the  smooth  sphere (primarily  specular  point  for  spheres 
with  large  kQD)  due  to  surface  roughness.  The  diffuse  scattering  contri¬ 
bution  to  the  phase  matrix  [S']  due  to  surface  roughness  is 

0  0 


[sD]  = 


[S?0  [s?J 
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12- 


[si]  [si] 
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[si]  [si] 


33 

[sj3]  [si,] 


’3k- 

.D  - 

sw»- 


(10) 


where 


[£”j] '  p[<0‘3V 


(11) 


in  which  A  =  itD  A  is  the  cross  sectional  area  of  the  unperturbed 

y 

particle  and  in  (ll)  is  the  like  (i=j)  or  cross  (i^j)  cross  polarized 

normalized  scattering  cross  section  (Bahar  and  Chakrabarti  1985).  The 
first  and  second  superscripts  i,J  denote  the  polarization  (V  vertical, 

K  horizontal)  of  the  scattered  and  incident  fields  respectively 

<olJ>  =  /  /  Jk  DlJ|2P0Q  sinydydfi/TT2  (12) 

D  J  1  o  2 

o  o 

in  which 

OO 

1 2 , 


Q  =  /  (x2(v»ar)  -  |x(v*ar)|  )exp(iv*rd)dxddzd 


(13) 


The  surface  height  autocorrelation  function  <hh*>  and  therefore  the  Joint 

2  2  ^5 

characteristic  function  Xo  (9)  are  only  functions  of  distance  r,5(x> z,) 

2  add 

measured  along  the  surface  of  the  unperturbed  sphere  since  the  random 
rough  surface  is  assumed  to  be  isotropic  and  homogeneous.  Furthermore, 

[SiiJ  =  Ay  i  <°21>D];,/,,,rp[ot3  {lk) 

for  i=3  (upper  sign)  and  i=U  (lower  sign).  For  if*J 

[sjj]  =  Ay  p[lm[±  <Cr22>D  +  (15) 


13 


■  .  *  h  V*  *  •  "  *1  *  -  *  *  •  w 


V-*  >  ,>  ^ 


Srfl 


vhere  the  upper  sign  is  used  for  i=^,  J=3  and  the  lower  sign  is  used  for 


i=3»  J-U.  In  (ll»)  and  (15) 

27T  7T 

<0k£>  =  /  /  kVJDki*P2Q  sinydYdfi/ir2  (16) 

"00 

In  (12)  and  (16)  P2  is  the  shadow  function  (Sancer  1969)-  It  is  the 
probability  that  a  point  on  the  surface  is  both  illuminated  and  visible 
given  the  slope  of  the  surface  at  that  point.  The  scattering  coefficients 
are  functions  of  n1,iif  and  ii  the  normal  to  the  unperturbed  surface 
of  the  particle  as  well  as  the  electromagnetic  parameters  of  the  material 
of  the  particle  e,y  (Bahar  1981).  Since  and  D1^  (i^j)  are  symmetric 
and  antisymmetric  respectively  with  respect  to  the  azimuth  angle, 6,  of  the 
particle,  the  remaining  eight  terms  of  [S^]  vanish. 

A  linearly  polarized  wave  is  assumed  to  be  normally  incident  upon  a 
parallel  layer  (of  optical  thickness  Tq),  containing  a  random  distribution 
of  irregular  shaped  particles.  Thus,  the  incident  Stokes  matrix  at 
z=0  (see  Fig.  l)  is 

<${|i'-l)6($ ' )  =  [l0]6(p,-l)6(<J>’)  (IT) 

in  which  6{*)  is  the  Dirac  delta  function.  Thus,  the  reduced  incident 
intensity  is 


IX  J  = 

inc 


A 


Clri]  ' 

The  (Uxl)  excitation  matrix  [1^3  in  (l)  is 

[Ij]  *  J  [S][Iri3dy *  [F]exp(-T ) 
in  which  [F]  can  be  expressed  as 
[F]  .  [S][I  ] 


y=i 

<J>'=0 


=  [F]  ♦  [F]  cos2$  +  [F]v sin24> 
0  a  0 


(18) 


(19) 


(20) 


and  the  matrix  [lo3  is  defined  in  (17).  The  solution  of  the  equation  of 


transfer  (l)  for  the  incoherent  specific  intensity  matrix  [i]  can  be 
expanded  in  a  Fourier  Series  (Cheung  and  Ishimaru  1982).  In  view  of  the 
excitation  (20),  only  the  azimuthally  independent  term  (m=o)  and  the 
c os2^  and  sin 2<p  terms  (m=2)  are  non-vanishing.  Thus,  the  transfer 
equation  (l)  together  with  the  associated  boundary  conditions 

[I]  =  0  ,  for  0  <  y  <  1  at  t  =  0  (2l) 

and 

[I]  =  0  ,  for  0  >  u  >  -1  at  T  *  T  (22) 

-  -  o 

are  solved  for  the  specific  incoherent  diffuse  scattering  intensities 
using  the  Gaussian  quadrature  method  (to  discretize  the  angles  6=cos  ^y) 
and  the  matrix  characteristic  value  technique  (Ishimaru  1978). 

The  only  remaining  parameter  that  needs  to  be  evaluated  to  solve  the 
equation  of  transfer  (l)  is  the  extinction  coefficient  (total  cross 
section)  The  "forward  scattering  theorem"  (Born  and  Wolf  196*0 

relates  the  total  cross  section  to  the  imaginary  part  of  the  scattering 

— r  -i 

(complex)  amplitude  in  the  forward  direction  f^(n  ,n  )  (Ishimaru  1978). 
Thus,  for  a  spherical  particle,  the  normalized  total  cross  section 
(per  unit  cross  sectional  area  of  the  sphere  of  radius  a)  is 

ot  =  — —  Im[f..(Hf,n1)]  i=l  ,2  (23) 

k  a 
o 

For  particles  with  irregular  shape  and  particles  with  random  rough 

surfaces,  it  is  very  difficult  to  determind  the  scattering  amplitude  f^. 

When  the  surface  roughness  of  the  particle  is  small,  such  that  the  roughness 
2  2 

parameter  £=lko<h  >  «  1,  perturbation  theory  can  be  used  to  determine  the 
diffuse  scattering  contribution  (Rice  1951,  Ruck  et  &1.  1970).  In  this 
case  (£  «  l)  0t  can  be  approximated  by  its  unperturbed  value  since  the 
diffuse  scattering  contribution  [S^]  for  large  particles  vanishes  in  the 
forward  direction  and  the  shadow  boundary  is  practically  the  same  as  that 
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of  the  unperturbed  particle.  However,  for  the  small  perturbation  case. 


the  effects  of  particle  surface  roughness  upon  the  specific  intensities 
is  almost  negligible.  Since  in  this  work  particles  with  very  rough 
surfaces  1  <  g  <  1*0  are  considered, small  perturbation  theory  cannot  be 
used  to  determine  the  scattering  amplitudes  f  and  the  effects  of  surface 
roughness  upon  the  specific  intensities  is  not  negligible.  To  facilitate 
the  solution  to  the  problem  when  B  is  large.  Judicious  use  is  made  of  the 
scattering  theorem  and  the  perceptive  observation  that  for  large  scatterers 
(physical  dimensions  large  compared  to  wavelength)  the  forward  scatter 
"shadow  forming  wave  is  the  same  for  all  surfaces  which  have  the  same 
shadow  line"  (Morse  and  Feshbach  1953).  Thus,  the  albedos  A  for  the 
particles  with  random  rough  surfaces  are  evaluated  as  follows 
o  a  a  a  ro.  \ 

s  s  to  .  s  to  I.  \ 


o_  -  s  toj 
t  ato  ■  atn 


a.  a 

to 


In  (2k)  a^,  the  normalized  scattering  cross  sections  for  the  finitely 
conducting  particles ,  are  given  by 


%  -  "T  /  lx2|°Kdn  *  IT  l<aV 


5  °sl  +  °s2  l25) 

in  which  X  is  the  rough  surface  characteristic  function  au  is  the  Mie 

solution  for  the  smooth  (unperturbed)  sphere  (Ruck  et  al.  1970,  Ishimaru 
1978)  and  dii  is  the  elementary  solid  angle.  Thus  integrate  with  respect 
to  the  azimuth  angle  to  get  (in  the  notation  used  by  Ishimaru,  1978) 

°sl  =  (V^/lxftlVe)!2*  |S2(e)|2]sin6d6  (26) 

The  full  wave  approach  is  used  to  obtain  the  diffuse  contribution  a g2  to 
the  scattering  cross  section.  Thus 

os2=0.25/(<aW>D-*-<oVH>D+<oHV>:D+<oHB>D)sinede  (27) 

in  which  are  the  like  and  cross  polarized  diffuse  differential 


/  ■*  .•  V  -  v , 


a  *  «.*•-»  ■«  a  > 

•a* 

>  *.*  'a  V- 


, ^  m'  -  .V  k'*  a'*  , 

■  *•  •  -  -  -  v  V*  -  »•  aViNJV-  ".Nk'j 


scattering  cross  sections  per  unit  solid  angle  (Bahar  and  Fitzvater  1983, 


Bahar  and  Chakrabarti  1985)-  Furthermore,  in  (2U),  o.  is  the  total 

to 

cross  section  for  the  smooth  particle.  On  approximating  (a^Q/a  )  by 
its  corresponding,  perfectly  conducting  value  (a  /o  )_  n  use  ma^e 

tO  X  Jr  •  v  • 

of  the  forward  scattering  theorem  and  the  fact  that  for  k^a  »  1,  the 
forward  scatter,  shadow-forming  wave  is  the  same  for  all  surfaces  which 
have  the  same  shadow  line.  Since  for  perfectly  conducting  particles 

^°t^P.C.  =  (asV.C.  ^28' 

where  (a  )  is  the  normalized  scattering  cross  section  for  the 

S  ±  •  0  * 

perfectly  conducting  particle  with  the  same  rough  surface  as  the  one 

under  consideration.  Thus  (a  )_.  n  is  evaluated  using  (25)  for  the 

s  r.L . 

corresponding  perfectly  conducting  particle  and  the  expression  for  the 
albedo  is  given  by 

A  ’  C<’Ao3/IlV®to3F.C.  *  VA2  •  (29: 

Two  limiting  cases  for  the  above  expression  are  of  particular  interest. 
As  the  conductivity  of  the  particle  increases  to  infinity  A^  -*  A^  and 
A  •*  1.  As  the  surface  roughness  becomes  very  small  B  -*•  0 ,  A^  -*•  1,  (since 


^°sV.C.  "*  ^°tV.C. 


(cr  )„  _ )  and  A  A,  *0  /(T  .  Thus  as  expected 
to  i  • U •  X  s  to 


A  -*•  1  for  perfectly  conducting  particles  and  A  -*•  A^  its  corresponding 

unperturbed  value  for  S'*  0.  The  value  of  the  extinction  coefficient 

is  obtained  from  the  relationship  0  =  O  /A. 

X  s 


p .. .. 


' .V.-  ; .  .... . . ... ... .. ... , 


4.1  Illustrative  Examples 

For  the  illustrative  examples  presented  in  this  section,  the  random 

rough  surface  height  h  (measured  normal  to  the  surface  of  the  unperturbed 

spherical  surface  of  the  particle  of  radius  a)  is  characterized  by  its 

surface  height  spectral  density  function  W(k  ,k  ).  Since  the  random  rough 

X  z 

surface  is  assumed  to  be  isotropic  and  homogeneous,  the  spectral  density 

2  2  k 

function  is  expressed  as  a  function  of  1^=  (k  +k 

A  X  z 

oo 

W(kT)  =  W(k  ,k  )  =  /  <hh,>exp(ikxxd+ik2zd)dx(jdzd 


=  -  /  <hh’>J  (k  r  )r  dr 
7T  '  o  T  d  d  d 

o 


In  (29)  J^k^r^  t^ie  Bessel  function  of  the  first  kind  and  <hh’>  is 

the  rough  surface  height  autocorrelation  function  which  is  only  a  function 

2  2 

of  distance  rd  =  (xd+zdPmeasured  along  the  surface  of  the  unperturbed 
particle.  In  view  of  the  Fourier  transform  relationship  between  W(k^)  and 


<hh’> 


<hh’>  =  / 


W(k  ,k  ) 


exp(-ik  x.-ik  z,)dk  dk 
x  d  z  d  x  z 


- 1  /  "WVd’W 
0 


The  following  explicit  expressions  are  assumed  in  this  work  for  the 


surface  height  spectral  density  function 


M  V  it 


_  r  k_  in 

C  T 

IT  -2.2 
k_+k 
_  T  m 


,  kT  >  0 


in  which  the  exponent  is  n  =  6  case  (a)  and  n  =  8  case  (b).  The  value  of 

W(k_)  is  maximum  for  km  =  k  (the  dominant  roughness  scale) 

I  i  m 

w(k  )  =  2C/TT  kn  (32) 

m  m 

For  n  s  8  the  scales  of  roughness  are  more  tightly  bound  around  k  than 

ZD 

2  2 

for  n  =  6.  The  mean  square  height  <h  >  and  mean  square  slope  <a >  of  the 


<h2>  =  I  /  w(kT)kTdkT  = 


and 


w 

<of>  =|  /  w(kT)k|dkT  = 


and  the  constant  C.  Thus 

c 

C/210  k 

JJJ 

,  n  *  8 

h 

CAOk 

m 

,  n  =  6 

(33 

C/81*  k1* 

m 

,  n  =  8 

• 

C/10  k2 
m 

,  n  =  6 

(31* 

and  the  surface  height  correlation  distance  r  (where  <hh'>  =  exp(-l)  is 


r  =  2(<h  2>/<02>)?5  =  ■ 


1.265/k  ,  n  *  8 

m 


l/k, 


m 


(35 


The  autocorrelation  functions  corresponding  to  the  rough  surface  height 
spectral  density  functions  (3l)  can  he  written  in  closed  form.  Thus  for 
n  =  8 


2  h  J>  "1 

R(0  =  [l  -  +  ^§2  +  307||CKl(O 


2  I»" 

k  96 


C2kq(0 


(36c 


ar.d  for  n  ®  6 

R(S 


,  _  f  3c2  nrK 
;  -  L1  -  ~ir  -  9srKi 


(O 


fi  +  4 
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C2ko(5) 


(36h 


in  which  K  and  Kn  are  the  modified  Bessel  functions  of  the  second  kind 
o  1 


or  order  zero  and  one  respectively  (Abramowitz  and  Stegun  196k)  and  the 
dimensionless  argument  C  is 

5  =  k  r.  (37; 

m  a 

The  wavelength  of  the  optical  excitation  is 

X  =  0.555m  (38! 

and  the  diameter  of  the  unperturbed  sphere  is 

D  =  10X  (39; 

The  particles  are  assumed  to  be  made  of  (i)  aluminum  for  which  the 


relative  dielectric  coefficient  is  (Ehrenreicb  1965) 


and  (ii)  a  dissipative  dielectric  vith 


(l*Ob) 


e  *  1.5-iB 
r 

Two  sets  of  data  for  each  case  n  =  6  and  n  *  8  and  each  dielectric 

material  are  presented  graphically  in  this  work.  In  the  first  set 

k  D  *  1*;  thus  the  correlation  length  r  is  fixed  (35)  while  both  the 
m  c 

2 

mean  square  height  and  slope  vary.  The  maximum  value  of  <og>  is  0.101 

corresponding  to  8ffi  =  l*k^<h^>  =  ^0  for  n  *  8  and  -  25  for  n  =  6 . 

2 

In  the  second  set  the  mean  square  slope  <Og>  is  set  equal  to  0.101  while 
both  the  mean  square  height  <h  >  and  k^  (and  thus  rc)  vary. 

In  Tables  T  through  VI  the  value  of  the  scattering  cross  sections  o 
a  _  and  c  (2k)  as  well  as  A,  ,A_  the  albedo  A  (28)  and  the  extinction 

S2  S  12 

coefficient  a  are  listed  for  different  values  of  the  roughness  parameter 

v 

2  2 

8  =  l*h  <h  >.  In  Tables  I  and  II  the  relative  dielectric  coefficient  is 
o 

£  =  -U0-il2  and  for  Tables  III  through  VI  t  -  1.5-i8.  Tables  I  through 

r  r 

IV  are  for  n  =  8  while  V  and  VI  are  for  n  =  6.  In  Tables  I,  III  and  V 

2  2 

k  D  =  4  and  <o  >  varies  while  in  Tables  II,  IV  and  VI  <o  >  =  0.101  and 
ms  s 

k  (and  the  correlation  distance  r  )  varies, 
m  c 

In  Fig.  2  the  scattering  cross  sections  og ,  Og^  and  o g^  are  plotted 

as  functions  of  8  (0  <  8  <  8  )  and  in  Fig.  3  the  albedo  A  and  the 

m 

extinction  cross  section  o  are  plotted  as  functions  of  8  for  the  case 

w 

given  in  Table  I.  Similarly,  Figs,  k  and  5  correspond  to  the  case  given 
in  Table  II,  Figs.  6  and  7  correspond  to  the  case  given  in  Table  III  and 
Figs.  8  and  9  correspond  to  the  case  given  in  Table  IV. 

It  is  interesting  to  note  that  while  Og  the  scattering  cross  section 
is  primarily  dependent  on  the  roughness  parameter  6,  the  absorption  cross 
section  o  =0.-0  depends  both  upon  8  and  the  mean  square  slope 

A  X  S 

2 

<0g>.  Thus  for  particles  with  the  same  roughness  parameter  6,  the 


.  «.*.  »• 


absorption  cross  section  a  increases  (and  the  albedo  decreases)  as 

a 

the  near,  square  slope  increases.  In  the  next  set  of  2h  plots,  the 
(differential)  like  and  cross  polarized  normalized  diffuse  scattering  cross 
sections  <o^>^  are  plotted  as  functions  of  the  angle  6  between  the 
incident  and  scattered  wave  normals  n1  and  n  respectively.  In  Figures  10 
through  18  the  s-urface  height  spectral  density  function  W  is  given  by 
(31)  vithn  =  8.  The  constant  parameter  for  the  different  curves  is 

p  2 

B  =  <h  >.  In  this  set  of  curves  k  D  =  U  thus  the  correlation  length 
o  m 

rc  is  constant  while  the  mean  square  slope  <a2>  is  different  for  different 

values  of  8.  The  relative  dielectric  coefficient,  for  Figures  10,  11  and 

12  is  e  =  -10-112  for  Figures  13,  1^  and  15  it  is  E  *  1.5-i8,  and  for 
r  r 

Figures  l6,  17  and  18  the  particle  is  assumed  to  be  perfectly  conducting 
( I  c  I  -*■  °°) .  In  Figures  19  through  27  n  =  8  however  the  mean  square  slope 

p 

is  constant,  <a  >  =  0.101  thus  for  different  values  of  the  roughness 
s 

raraneter  B  the  correlation  distance  r  is  different.  In  Figures  28 

c 

PG 

through  33  <(?  is  given  for  n  =  6  and  e  =  1.5-i8.  For  Figures  28, 

29  and  30,  k  D  =  1  while  the  mean  square  height  and  mean  square  slope 
m 

2 

are  parameters  and  for  Figures  31,  32  and  33  <0  >  =  0.101  while  the  mean 

square  height  and  correlation  distance  rc  are  parameters. 

In  order  to  facilitate  comparison  between  the  different  data  for  the 

case  n  =  8,  the  data  corresponding  to  B  =  ^0  is  common  to  the  set  k  D  =  ^ 

m 

2 

and  the  set  <ag>  =  0.101.  Similarly,  for  the  case  n  =  6,  the  data 

corresponding  to  B  =  25  is  common  to  the  set  km  D  *  U  and  the  set 

<a2>  *  0.101. 

s 

2 

For  the  case  <0  >  constant  the  diffuse  scattering  terms  merge  for 
s 

B  >  5  except  for  scattering  in  the  near  forward  direction  where  the 
parameter  B  (and  r )  has  a  very  significant  effect  on  <o2^>^.  For  the 
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2 

case  k  D  =  ^  (r  =  const),  the  parameter  8  (and  <c  >)  has  an  effect  on 
me  s 

po 

the  diffuse  scattering  cross  sdetion  <0  for  all  6. The  main  differences 
are  in  the  backscattered  (8=0)  and  the  forward  scattered  (6=l80°) 
directions.  The  effects  of  particle  conductivity  (absorption)  on  the 
scattering  cross  sections  are  also  evident  from  the  plots  in  these  figures. 
The  level  of  the  diffuse  scattering  terms  decreases  as  the  absorption 
of  the  particle  increases.  The  differences  between  the  results  for 
cases  n  =  8  and  n  =  6  are  most  significant  for  scattering  in  the  forward 
direction. 

Note  in  particular  that  there  is  a  significant  like  polarized  back- 

scatter  enhancement  due  to  the  surface  roughness  of  the  particle.  The 

cross  polarized  diffuse  cross  sections  have  sharp  nulls  in  the  backscatter 

direction.  This  is  due  to  the  symmetry  of  the  particle  and  the  fact  that 

the  particle  surface  roughness  is  assumed  to  be  isotropic. 

In  the  last  set  of  figures,  solutions  of  the  transfer  equation  (l) 

are  presented  for  a  parallel  layer  consisting  of  a  random  distribution  of 

particles  with  rough  surfaces.  The  incident  wave  is  linearly  polarized 

and  normally  incident  upon  the  layer  of  optical  thickness  T  .  The  special 

cases  considered  in  detail  _  -e  for  particles  with  very  rough  surfaces 

2 

8  =  1*0,  k  D  =  1*,  <G  >  =  0.101,  n  =  8,  X  -  0.555y-  Both  aluminum  particles 
m  s 

=  — h 0— i 12  and  dissipative  dielectric  particles  =  1.5-i8  are  considered 
In  Figures  3^  and  35  the  transmitted  diffuse  incoherent  specific 
intensities  I  (copolarized)  and  Ig  (cross  polarized)  are  plotted  as 
functions  of  0  (0  <  6  <  90°)  and  in  Figures  36  and  37  the  reflected  diffuse 
specific  intensities  are  plotted  as  functions  of  0  (90°  <  0  <  l80°).  The 
optical  thickness  of  the  layer  is  *  2  and  the  relative  dielectric 
coefficient  is  =  -l<0-il2.  Single  scatter  and  multiple  scatter  results 
are  shown  for  both  (except  for  single  scatter  cross  polarized  intensity 
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for  the  smooth  particles).  Similarly,  in  Figures  38  through  1*1,  the 
copolarized  and  cross  polarized  transmitted  and  reflected  specific 


intensities  are  plotted  as  functions  of  6  for  T  =  2  and  e  =  1.5-i8. 

or 

In  Figures  1*2  through  1*7  the  copolarized  incoherent  specific  intensity 


1^  (scattered  electric  field  parallel  to  incident  electric  field)  and  cross 


polarized  incoherent  specific  intensity  I  (scattered  electric  field 

y 


perpendicular  to  incident  electric  field)  (ishimaru  1978)  are  plotted  as 


functions  of  the  azimuth  angle  4>-  The  optical  thickness  Tq  =  2  and 


a  -I4O-H2.  In  Figures  1*2  and  1*3,  1^  and  I  are  plotted  for  0  =  9*7  - 


Similarly  in  Figures  1*1*  and  1*5,  I  and  I  are  plotted  for  0  =  170.3  ,  and 

x  y 


in  Figures  1*6  and  1*7  I  and  I  are  plotted  for  0  =  175-8  .  Here,  too,  both 

x.  y 


the  single  scatter  and  multiple  scatter  results  are  given  for  both  the 


smooth  and  the  rough  particles. 

In  Figures  1*8  through  50,  the  degree  of  polarization  m  (Ishimaru  1978) 

is  plotted  as  a  function  of  the  azimuth  angle  <J>.  The  optical  thickness 

is  T  =2  and  the  relative  dielectric  coefficient  is  e  =  -l*0-il2.  In 
o  r 

Figure  1*8  0  =  l*. 2°  and  9-7°,  in  Figure  1*9  9  *  15-3°  and  20.8°  and  in 
Figure  50  6  =  170.3°  and  6  =  175-8°. 

1*.2  Computer  Data  and  Programs 

The  tapes  contain  all  the  computer  programs  needed  in  this  work  as  well 

as  the  computed  data  for  the  following  12  cases.  Surface  height  spectral 

density  function  W  (31)  with  n=6  and  n=8.  Complex  dielectric  coefficient 

ey  =  -l*0-il2,  er=  1.5-i8  and  (e^J  -*■«*>.  The  mean  square  slope  <o^>  is  a 

constant  parameter  (0  and  r  vary)  in  one  set  of  data  while  the  correlation 

c  2 

length  r  is  a  constant  parameter  (0  and  <0  >  vary)  in  the  other  set  of 
c  s 

data.  Data  for  layers  of  different  optical  thickness  =  0.1,1,2,10 
are  computed.  Selected  data  are  plotted  in  Figures  2  through  50. 
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TABLE  II.  Extinction  Cross  Section  and  Albedos 
for  Spheres  vith  Rough  Surfaces 
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8.  Figure  Captions 
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Fig.  1.  Scattering  geometry  indicating  incident  and  scatter  vave  normals 


Fig.  2. 

Fig.  3. 

Fig.  1*. 

Fig.  5. 

Fig.  6. 

Fig.  7. 

Fig.  8. 


n1  and  n  and  corresponding  field  components  E^  parallel 
(vertical)  and  perpendicular  (horizontal)  polarizations. 

Scattering  cross  sections  O  0g2  and  0g  versus  the  roughness 
parameter  8  (Table  l). 

"Extinction  cross  sections  0  and  albedo  A  versus  the  roughness 

v 

parameter  B  (Table  I). 

Scattering  cross  sections  0  ,  ,  0  _  and  0  versus  the  roughness 

Sl  S^:  S 

parameter  B  (Table  II ). 

Extinction  cross  sections  0^  and  albedo  A  versus  the  roughness 
parameter  B  (Table  II). 

Scattering  cross  sections  0  . ,  0  _  and  o  versus  the  roughness 

si.  s^.  s 

parameter  8  (Table  III). 

Extinction  cross  sections  0  and  albedo  A  versus  the  roughness 
parameter  8  (Table  III). 

Scattering  cross  sections  c  , ,  0  „  and  o  versus  the  roughness 

si  s2  ° 


parameter  8  (Table  IV). 


Fig.  9-  Extinction  cross  sections  0^  and  albedo  A  versus  the  roughness 
parameter  8  (Table  IV). 
w 

Fig.  10.  <0  >  versus  0,n=8,  k  D  =  1* ,  e  =  -l*0-il2 

°  D  m  r 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  ho  (0). 

HH 

Fig.  11.  <0  >  versus  0,  n  =  8,  kD  =  l*,  e  =  -l*0-i!2 

°  D  m  r 

B  =  1  (0),  5  (A),  10  (+),  20  (x),  hO  {O). 

VK  KV 

Fig.  12.  <0  >£  =  <0  >D  versus  0,  n  =  8,  k^D  =  Cr=  -l*0-il2 

8=1  (0),  5  (A),  10  (  +  ),  20  (x),  hO  (O). 
w 

Fig.  13.  <o  >D  versus  0,  n  =  8,  k^D  =  1*,  Cr=  1.5-i8 

B  =  1  (0),  5  (A),  10  {+),  20  (x),  40  (O). 

HH 

Fig.  ll*.  <0  >p  versus  0,  n  =  8,  k^D  =  1*,  er=  1.5-18 
B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

VH  HV 

Fig.  15.  <0  >n  =  <0  versus  0,  n  =  8,  k  D  =  1*,  e  =  1.5-i8 

VO  m  r 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  «». 
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Fig.  16.  <0  >D  versus  Q,  n  =  8,  kmD  =  1*,  perfect  conductor 

B  =  1(0),  5  (A),  10  (+),  20  (x),  1*0  «>)• 

KK 

Fig.  17.  <0>D  versus  0,  n  =  8,  k^D  =  1* ,  perfect  conductor 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

VH  TTV 

Fig.  18.  <0  >„  =  <0  >  versus  6,  n  =  8,  k  D  O,  perfect  conductor 

D  D  ID 

6  =  1  (0),  5'  (A),  10  (  +  ),  20  (x),  1*0  (£>). 

W  2  i 

Fig.  19.  <o  >D  versus  0,  n  =  8,  <0g>  =  0.101,  £r=  -l*0-il2 

B  *  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  <0>- 

mi  p 

Fig.  20.  <0  versus  0,  n  =  8,  <0g>  =  0.101,  er=  -l*0-il2 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (<>)• 

Fig.  21.  <0VH>d  =  <oHV>d  versus  0,  n  =  8,  <0^>  =  0.101,  e,=  -l*0-il2 
B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

Fig.  22.  <o^>  versus  0,  n  =  8,  <0^>  =  0.101,  er=  1.5-i8 

B  =  1  (0),  5  (A),  10  (+),  20  (x),  1*0  (O). 

Fig.  23.  <0^*‘>D  versus  0,  n  =  8,  <0^>  =  0.101,  er=  1 . 5— i 8 

B  *  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (£>). 

Fig.  2l*.  <a^>p  =  <a^V>p  versus  6,  n  =  8,  <o^>  =  0.101,  er=  1.5-i8 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

W  „  2 

Fig.  25.  <0  >D  versus  0,  n  =  8,  <0g>  =  0.101,  perfect  conductor 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  <»• 

HK  2 

Fig.  26.  <0  >D  versus  0,  n  =  8,  <0g>  =  0.101,  perfect  conductor 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (<». 

VH  HV  2 

Fig.  27.  <0  >D  =  <0  >D  versus  0,  n  =  8,  <Os>  =  0.101,  perfect  conductor 

B  *  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

Fig.  28.  <c^>  versus  0,  n  =  6,  k  D  =  1* ,  e  =  1.5-i8 

D  ffir 

B  =1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (<^). 

ini 

Fig.  29.  <0  >p  versus  0,  n  =  6,  k^D  =  1*,  1.5-i8 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  «>). 

Fig.  30.  <0*%  =  <aHV>j)  versus  0,  n  *  6,  k^D  =  1*,  cT-  1.5-i8 

B  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

W  2 

Fig.  31.  <0  >D  versus  0,  n  =  6,  <0g>  -  0.1,  Ej.*  1.5-18 

B  «  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O)- 

*3  A 
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FF  2 

Fig.  32.  <o  ‘  versus  0,  n  =  6,  <Os>  =  0.1,  er=  1 . 5-i 8 
8  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (O). 

Fig.  33.  <0  >D  =  <c“  >p  versus  6,  n  =  6,  <a^>  =  0.1,  e^=  1 . 5— i  8 

8  =  1  (0),  5  (A),  10  (  +  ),  20  (x),  1*0  (C). 

Fig.  31*.  Specific  incoherent  intensity  I  for  a  linearly  polarized 

o 

Normal  incidence,  8  =  1*0,  kmD  =  1*,  <Og>  =  0.101,  n  =  8  ( 31 ) 

A  =  0.555V,  e  =  -l*0-il2,  D  =  10A.  Transmitted  T  =  2,  d>  =  0. 

r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (  +  )  smooth  (A)  rough. 

Fig.  35.  Specific  incoherent  intensity  (I  )  for  a  linear V  polarized  vave. 

^  p 

Normal  incidence,  8  =  1*0,  k^D  =  1*,  <a^>  =  0.101,  n  =  8  (31) 

A  =  0.555V,  £r  =  -l*Q-il2,  D  =  10A.  Transmitted  Tq  =  2,  <J>  =  0. 
First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (+)  smooth  (A)  rough. 

Fig.  36.  Specific  incoherent  intensity  1^  for  a  linearly  polarized  wave. 

Normal  incidence,  8  =  1*0,  k^D  =  1*,  <o^>  =  0.101,  n  =  8  (31) 

A  =  0.555v,  e  =  -l*0-il2,  D  =  10A.  Reflected  T  =  2,  4  =  0. 
r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (■*■)  smooth  (A)  rough. 

Fig.  37.  Specific  incoherent  intensity  (I  )  for  a  linearly  polarized  wave. 

^  2 

Normal  incidence,  8  =  1*0,  k  D  =  h ,  <0  >  =  0.101,  n  =  8  (3l) 

m  s 

A  =  0.555v,  e  =  -l*0-i!2 ,  D  =  10A.  Reflected  T  =2,<{)  =  0. 
r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (+)  smooth  (A)  rough. 

Fig.  38.  Specific  incoherent  intensity  I  for  a  linearly  polarized  wave. 

^  2 

Normal  incidence,  8  =  1*0,  k  D  =  1*,  <c  >  =  0.101,  n  =  8  ( 31 ) 

ms 

A  =  0.555V,  E  =  1.5-i8,  D  =  10A.  Transmitted  T  =  2,  4  =  0. 
r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (+)  smooth  (A)  rough. 

Fig.  39-  Specific  incoherent  intensity  1.  for  a  linearly  polarized  wave. 

^  2 

Normal  incidence,  8  =  1*0,  kffiD  =  1*,  <og>  =  0.101,  n  =  8  (31) 

A  =  0.555V,  c  =  1.5-i8,  D  =  10A.  Transmitted  T  =  2,  <b  -  0. 
r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (+)  smooth  (A)  rough. 
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1+0.  Specific  incoherent  intensity  for  a  linearly  polarized  wave. 

Normal  incidence,  8  =  1+0,  k  D  =  1+ ,  <a^>  =  0.101,  n  =  8  (3l) 

m  s 

A  =  0.555V,  er  =  1.5-i8,  D  =  10X.  Reflected  tq  =  2,  <J>  =  0. 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (+)  smooth  (A)  rough. 

1+1.  Specific  incoherent  intensity  (l5)  f°r  a  linearly  polarized  wave 

<-  o 

Normal  incidence,  8  =  l+0,kD  =  l+,  <0>  =  0.101,  n  =  8  (31) 

m  s 

X  =  0.555V,  er  =  1.5-i8,  D  =  10  .  Reflected  tq  =  2,  <t>  =  0. 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (+)  smooth  (A)  rough. 

1+2.  Specific  incoherent  intensity  I  for  a  linearly  polarized  wave. 
Normal  incidence,  8  =  1*0,  k  D  =  1+,  <0^>  =  0.101,  n  =  8  (31), 

IQ  S 

\  =  0.555V,  £-r  =  -^0-il2,  D=10X.  Transmitted  =  2 ,  6  =  9  -7  . 
First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (X)  smooth  (X)  rough. 

1+3.  Specific  incoherent  intensity  1^  for  a  linearly  polarized  wave. 

Normal  incidence,  8  =  1+0,  kD  =  l+,  <o^>  =  0.101,  n  =  8  (31), 

ms 

X  =  0.555V,  £  =  -l*0-il2,  D  =  10X.  Transmitted  T  =  2,  0  =  9-7  . 

r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (X)  smooth  (X)  rough. 

11- .  Specific  incoherent  intensity  I  for  a  linearly  polarized  wave. 

x  2 

Normal  incidence,  8=l+0,kD  =  l+,  <a>  =  0.101,  n  =  8  (31), 

m  s 

X=0. 0555V,  £  =  -1+0-H2 ,  D  =  10X.  Reflected  T  =  2,  6  =  170.3° 

r  o 

First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (X)  smooth  (X)  rough. 

1+5-  Specific  incoherent  intensity  1^  for  a  linearly  polarized  wave. 

Normal  incidence,  8  =  1+0,  k  D  =  1+,  <o^>  =  0.101,  n  =  8  (31), 

m  s 

X  =  0.555’J,  er  =  -l+0-il2,  D  -  10X.  Reflected  tq  =  2,  0  =  170.3' 
First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  {X)  smooth  (X)  rough. 


1+6.  Specific  incoherent  intensity  I  for  a  linearly  polarized  wave. 

*  2 

Normal  incidence,  8  =  l+0,kD  =  l*,<o>  =  0.101,  n  =  8  (31) 

m  s 

X  =  0.555V,  £r  =  -l+0-il2 ,  D  =  10X.  Reflected  xo  =  2,  6  =  175.8' 
First  order  ( — ),  smooth  and  rough  particles.  Multiple 
scatter  (X)  smooth  (X)  rough. 
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Fig.  J*7.  Specific  incoherent  intensity  I  for  a  linearly  polarized  wave. 

y  2 

Mormal  incidence,  B  =  1*0,  k  D  =  1*,  <o  >  =  0.101,  n  =  8  (31), 

ms 

X  =  0.555y,  e  =  -U0-112,  D  =  10X.  Reflected  T  =  2,  0  =  175. 
r  o 

First  order  ( — ) ,  smooth  and  rough  particles .  Multiple 
scatter  (X)  smooth  (x)  rough. 

Fig.  1*8.  Degree  of  polarization  m.  Normal  incidence,  linearly  polarized 

wave,  6  =  40,  k  D  =  1*,  <a2>  =  0.101,  n  =  8  (31),  X  =  0.555U, 
m  s 

Ers-k0-il2 >D  =  10X ,  =  2.  8  =  1*.2°  (+)  smooth  (X)  rough; 

0  =  9*7°  (0)  smooth  (A)  rough. 

Fig.  1*9.  Degree  of  polarization  m.  Normal  incidence,  linearly  polarized 

wave,  6  *  40,  k^D  =  1*,  <o2>  =  0.101,  n  =  8  (31),  X  =  0.555y, 

e  =-l*0-il2jD  =  10X ,  t  -2.  0  =  15*3°  (  +  )  smooth  (X)  rough 
r  o 

8  =  20.8  (0)  smooth  (A)  rough. 

Fig.  50.  Degree  of  polarization  m.  Normal  incidence,  linearly  polarized 
wave,  8  =  40,  k  D  =  4,  <o2>  =  0.101,  n  =  8  (3l),  X  =  0.555y, 

HI  S 

er ”l*0-il2jD  =  10A,  T  =2.  0  -  170.3°  (  +  )  smooth  (X)  rough 
0  =  175.8°  (0)  smooth  (A)  rough. 


APPENDIX  B 


0.95 


it' 


Fig.  II 

appendix  b 


o • uu  -Id. 00 


Fig.  20 


APPENDIX  B 


61 


Fig.  21 


PENDIX  B 


’  •  h  ^  ^  »  *  *  •  '  •  *  '  .  *  •  *'*.“.  *«•*.**■*  i.^  . 

'  %  *«  *«.  .  %  V  ..  A  /.  •*  «  '  -Oa*.  O.-J  *  *  •.*  A*  -  "  •  ' •  *  *  v*  l*  -  •  v'  »*  O.  •_-  . 


e 


Fig.  24 


APPENDIX  B 


65 


APPENDIX  B 


/  r* 


^-4 


3 


Fig.  32 


NDIX  B 


'  ' -v*  vVy-1  % 


v*\  ••  .  .  V 


smooth 


[»!•] 


iWiTi 


nmu 


9 


Fig.  38 


9 


appendix  b 


APPENDIX  B 


85 


